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Abstract 

We show a relation between fractional calculus and fractals, based 
only on physical and geometrical considerations. The link has been 
found in the physical origins of the power-laws, ruling the evolution of 
many natural phenomena, whose long memory and hereditary prop¬ 
erties are mathematically modelled by differential operators of non 
integer order. Dealing with the relevant example of a viscous fluid 
seeping through a fractal shaped porous medium, we show that, once 
a physical phenomenon or process takes place on an underlying fractal 
geometry, then a power-law naturally comes up in ruling its evolution, 
whose order is related to the anomalous dimension of such geometry, as 
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well as to the model used to describe the physics involved. By lineariz¬ 
ing the non linear dependence of the response of the system at hand 
to a proper forcing action then, exploiting the Boltzmann superpo¬ 
sition principle, a fractional differential equation is found, describing 
the dynamics of the system itself. The order of such equation is again 
related to the anomalous dimension of the underlying geometry. 


1 Introduction 

It is widely believed that the origin of the concept of fractional deriva¬ 
tive dates back to the 1695, when Leibniz, in its correspondence with de 
L’Hospital, posed the question about the meaning of a derivative of order one 
half. That intuition represented the seed for mathematicians of the following 
three centuries as Liouville, Riemann, Griinwald, Letnikov et ah, to develop 
the theory of the fractional calculus as it is nowadays known. However, de¬ 
spite of its long time origins, only in the recent few decades the fractional 
calculus gained the attention of the scientihc community, having been con¬ 
sidered, for most of the time, as a pure theoretical held of research. The rise 
of the interest on the fractional calculus, has to ascribe to its suitability in 
modelling a large variety of phenomena, ranging almost all the helds of the 
naturals sciences. Nowadays, its applications can be found in biology and 
biomechanics [1-6], in modelling the rheological properties of viscoelastic ma¬ 
terials [7-9], in electrical circuits [10], in anomalous transport and diffusion 
processes in complex media [11-14], and in many other branches of physics 
and engineering [15-17]. The powerful of the fractional operators relies on 
their non-local power-law structure, that makes them suitable in describing 
the memory and hereditary properties of many physical phenomena and pro¬ 
cesses, and in particular in modelling the “anomalous” dynamics of complex 
systems, where scaling power-laws emerge as a consequence of such complex¬ 
ity. It has been shown in fact, that fractional (integro-)differential equations 
naturally arise once systems characterized by power law long range inter¬ 
actions, or by power law memory, are considered [18-20]. Despite of their 
widely recognized usefulness in describing nature, never a clear and uni¬ 
versally accepted geometrical interpretations has been pointed out for such 
operators. 

The development of the theory of fractals, formulated by Mandelbrot in 
the second half of the twentieth century [21], who consolidated the work of 


2 



many other prior mathematicians and scientists in general, opened a new 
perspective to this aim. Since that, several attempts have been made, in or¬ 
der to elucidate a connection between the fractional operators and this new 
geometry, whose anomalous non-integer dimension seemed to be the key of 
the story. The concept of fractance has been one among these attempts. 
By fractance it is usually referred in literature, a mechanical (or electrical) 
model, in which classical elements such as springs and dashpot in the case of 
viscoelasticity (or resistors and capacitors in electrical circuits) are disposed 
as a fractal tree, with an opportune power law scaling postulated for the rel¬ 
evant physical properties, in order to obtain a constitutive law of non-integer 
order for the whole system [4, 10, 22, 23]. Another line of research has been 
developed by Nigmatullin and Le Mehaute [24, 25], who showed that the 
procedure of averaging of a smooth function over a non dense self-similar 
set, leads to a Riemann-Liouville integral operator, whose order is related 
to the anomalous dimension of the set itself. Working in the time domain, 
they proposed the interpretation of their result as manifestation of a partial 
(fractal) memory of the system. A different approach, that led to the dehni- 
tion of a new kind of operators, acting on a non dense support and known 
in literature as local fractional operators, has been pursued by Kolwankar 
and Gangal [26-32]. The idea underlying their work is that, as they are 
dehned on a continuum support, the ordinary fractional operators cannot 
be dehned on a non dense object like a fractal. This approach, however, 
required the dehnition of proper local operators, different from the classi¬ 
cal Riemann-Liouville (R-L) or Caputo’s fractional derivatives and integrals, 
missing then the aim of giving a connection between fractal geometry and or¬ 
dinary fractional calculus. Many others attempts have been proposed in the 
last decades. The property of fractional operators of changing the anomalous 
dimension of a random (or deterministic) fractal object, to which they are 
applied, has for example been shown in [33]. A similar result was also derived 
by Rocco and West in [34], where fractal functions have been shown to have 
fractional derivatives and integrals, whose anomalous dimension differs from 
the original one. Moreover, by following a different approach, the authors 
recovered the definition of a critical order of derivation and integration, de¬ 
duced by Kolwankar and Gangal in [30, 31]. However, in spite of such effort, 
never a clear relation between fractals and fractional calculus, only based on 
geometrical considerations, has been pointed out up to now. 

It is our belief, that such a link has not to be sought in the re-de£nition 
of fractional operators on fractal support, or in the construction of ad-hoc 
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mechanical (or electrical) models, but in the physical origins of the scaling 
power laws that emerge from the dynamics of a large variety of complex 
systems, that are at the basis of the long memory and hereditary properties 
commons to many natural phenomena, and that is the main peculiarity of 
fractional operators. In this paper, we want to show that such power laws 
come out in a natural way, when physical phenomena taking place on an 
underlying fractal structure are considered. It turns out that the order of 
such power laws is strictly related to the dimension of the geometry on which 
the phenomena take place, as well as to parameters related to the physics 
of the problem. The link between fractals and the raising of dimensional- 
related power laws is not new in literature and has been investigated in the 
past, with particular emphasis in diffusion phenomena in complex media 
[35-38]. It has been shown for example that, modelling the dynamics of 
diffusion processes in terms of a random walk, the time dependence of the 
mean square displacement away from the initial position, deviates from the 
linearity predicted by the Tick’s law for diffusion in homogeneous media, 
and a power-law of non integer order, related to the fractal dimension of the 
medium, has been found. 

The reason of the emergence of such power laws has to be sought in 
the dehnition of fractals as self-similar structures, that differs from regular 
geometric shapes in their fractional dimensional scaling. In other words, a 
(deterministic) fractal can be thought as an object, each part of which is a 
scaled copy of the whole, that differ from a regular geometrical shapes in its 
non integer scaling coefficient. If we rescale by a factor A the linear dimension 
of a square or a cube, for example, the corresponding area and volume, scale 
as A^ and A^, because D = 2 and D = 3 are their respective dimensions. On 
the other hand, if we rescale the linear dimension of a fractal object by the 
same factor, its spatial content scales by the factor A", where D = a is its 
anomalous dimension. As we show in the present paper, we may expect that, 
if a physical phenomenon or process takes place on an underlying self-similar 
structure, then all the scales have to manifest in its evolution, along with the 
relative scaling factor. In other words, the functional spatial and/or tem¬ 
poral dependence, ruling the phenomenon at hand, should manifest scaling 
properties related to the ones characterizing the underlying geometry. 

The connection between fractals and the emergence of non integer order 
power laws, appears immediately clear, once the self-similar structure of the 
power type functions is pointed out. If we zoom out by a factor A the time 
scale in the power law f(t) = t°‘, the rescaled function reads as: /(At) = A" t“. 
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that still is the same power law, rescaled by the factor A“. The self-similarity 
and scaling property of the are then exactly the same as for a fractal 
object of anomalous dimension D = a. From what said, we should expect 
that the order of the observed power laws, ruling the evolution of the natural 
phenomena, are strictly related to the dimension of the underlying geometry 
(as well as to the physics of the problem at hand). 

Because of such considerations, aim of the paper is, in first instance, 
to convince the reader (with the aid of relevant examples) that the previous 
arguments hold true for whichever phenomenon or process taking place on an 
underlying fractal geometry, because of the generality of the geometrical and 
physical considerations on which they are based. Moreover, the main goal of 
the paper is to perform a step forward and £11 the gap between fractals and 
fractional calculus. Basing on the relationship established between power- 
laws and fractal geometries, we will show that, once we assume sufficiently 
small the amplitude of the loads acting on the system, then a fractional 
differential equation naturally comes up, ruling the dynamics of the system 
itself, whose order turns out to be related to the dimension of the underlying 
geometry. 

In the following sections, we will show the validity of these considerations 
studying the time evolution of the flow of a viscous fluid through a porous 
medium. As expected, a power-law time dependence for the overall flux has 
been found, whose order is related to the anomalous dimension of the pattern 
along which we suppose the fluid flows, and to a coefficient related to the 
physical model by which the motion of the fluid itself is described. 

2 The main model 

In this and in the following sections, we prove the emergence of a power law 
time dependence in the flux of the viscous fluid flowing through the porous 
medium. We find an order for the power-law, related to the anomalous di¬ 
mension of the fractal shape by which we model the pattern of the streams 
along which the fluid flows. Although we don’t give any rigorous mathe¬ 
matical demonstration of the validity of this statement independently on the 
physical situation at hand, we are confident that, because of the generality 
of the geometrical and physical considerations on which the following argu¬ 
ments are based, the results derived in the paper hold true for whichever pro¬ 
cess evolving on an underlying fractal structure. Fractal-like structures are 


5 




Figure 1: Sketch of a Lichtenberg figure, used to model the pattern according 
to which the fluid flows through the porous medium 


widely observed in nature, and the shapes of lightening, trees, river’s estuary, 
ferns, etc. are examples. Since that, the “tree” pattern shown in fig.l, known 
in literature as Lichtenberg hgure [39, 40], from the name of the physicist 
who originally discovered it in the discharging of insulating material, is likely 
to describe the flow of a viscous fluid through a porous medium. Such pat¬ 
tern consists of a main stream, branching off in different secondary streams, 
with progressively smaller diameter, as we go downwards from z = 0. In the 
limit —)■ oo, i.e. carrying forward the branching process up to inhnity, the 
cross section of such pattern reveals a fractal stochastic shape, in which each 
stream is surrounded by secondary streams. For simplicity’s sake, in our 
calculations, we will replace such complex structure with an hollowed brick, 
in which all the streams are equally long, and distributed according to a de¬ 
terministic fractal shapes such as the square or triangular Sierpinsky carpets 
shown in fig.2. This simplifying hypothesis doesn’t affect the physical mean¬ 
ing and the generality of our arguments, because we have simply replaced the 
complex structure in hg.l, by a much simpler one, but still having the key 
scaling properties characteristic of a fractal shape. In what follows we start 
by setting up the general calculations leading to an expression for the time 
evolution of the out coming flow, based only on geometrical considerations, 
and so valid for whichever physical situation we are interested in. In the 
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Figure 2: The square and triangular Sierpinski carpets 

next sections, we will address the following specihc problems: i) in section 
3, we start by proving the natural appearance of a power-law in ruling the 
time evolution of the flux through the brick dehned above, that we take as a 
model for the pattern according to which the streams, along which the fluid 
flows, are distributed. To this aim, we study the phenomenon of percolation 
of the fluid, assuming it is subjected only to the gravitational and friction 
forces, ii) we show in section 4 that, if a time-dependent load acts on the 
system, that we might identify in a net pressure history applied to the brick, 
then, once the hypothesis of linear regime is assumed, a constitutive law of 
non integer order emerges in relating snch pressnre to the time evolntion of 
the flow throngh the brick. 

Let us start by dehning the notation that will be used through the pa¬ 
per. We indicate the nnmber of n—order streams and the total area of the 
corresponding cross-section by and S'„ = respectively, being y 

a form factor characteristic of the cross section, N the branching factor, and 
C = its characteristic length (in which Iq is the qnantity relative to the 
primary stream, and 1/e is the scaling factor of the fractal shape considered). 
We indicate in general with Vn{t) the velocity of the flnid, flowing along the 
n—th order stream, which depends on the physical model used to describe 
the motion of the flnid as well as on the scaling law of the fractal at hand. We 
will carry ont the calculations in the continnnm limit for which n —)■ a G M, 
that means to treat the discrete fractal as a continnnm one, so that the above 
dehned qnantities become fnnctions of the variable a. Pnshing the branching 
process np to inhnity, the ont coming how 0(t) of the hnid, howing throngh 
the discrete fractal shape, has the form: 

OO 

0(t) = y^^SnNnVnjt) (1) 

n=0 

in which the snm runs over all the orders of the fractal shape, according 
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to which the streams are distributed. The same equation, in the continuum 
limit, reads as: 

roo 

= / da S{a)N{a)v{a,t) (2) 

Jo 

Such equation can be rewritten in therms of only the (inverse of the) scaling 
factor e, and the anomalous dimension D = 7 ^ of the fractal at hand, as: 

’ me ’ 

roo 

= X^l / dae^^~'^'>°‘v{a,t) = x^o^efr(t) (3) 

Jo 

having dehned the effective velocity 

poo 

Ves{i) = / da (4) 

Jo 

Since the velocity of the flow depends on the order a only through the 
corresponding characteristic length i{a) = it is useful to perform the 

change of variable 

X = e“"; dx = —x In e da] a = Q ^ x = 1] a = +cxd —)■ a; = 0 (5) 

and so to write eq.(4) as 

= r-— [ dx x^^~^^v'{x,t) (6) 

log e Jo 

being v'{x,t) the functional form of the velocity, in the new variable x. 
The effective velocity is so given by the sum of all the stream’s velocities, 
weighted by the power-type kernel x^~^. It is worth to stress that such 
kernel depends only on the dimension D of the shape by which we describe 
the cross section of the fractal brick and, as it is evident from eq.(4), in the 
integer limit D —)■ 2 the effective velocity is given simply by the sum of the 
velocities of the different streams , because in this limit the cross section of 
the brick scales as a bi-dimensional object, and so the total area relative to 
the streams of whichever order, is equal to the area of the main one. 

Up to now, nothing has been said about the physics involved, that is 
encoded in v{x, t), but we obtained eqs.(4) and (6) only basing on geometrical 
considerations. In the next sections we will introduce physics in our model, 
describing the flow along each stream by using the Bernoulli equation, and 
modelling the losses by the Darcy-Weisbach law. As anticipated, we will 



show that in both the situations considered, a power law for the overall flux 
naturally conies up, whose order is related to the anomalous dimension of 
the pattern according to which the streams are distributed (i.e. to the fractal 
dimension of the cross section of the brick). 

3 Power laws from fractals: the percolation 
problem 

Let us start by considering the case in which the brick is saturated with the 
fluid, at rest for f < 0, and let us study the time evolution of the out coming 
flux for t > 0. The only forces that we suppose acting on the system in 
such conhguration, are the gravitational and the friction ones. Assuming the 
area of each stream small enough so that the friction forces are predominant 
respect to the gravitational ones, the velocity of the fluid can be considered, 
with good approximation, constant along each stream, so that we can write: 

v'{x, t) = V{x, h) [H{t) — H{t — T{x, h)] (7) 

where H{x) is the unit step function, and V{x, h) and T(a;, h) are respectively 
the constant value of the velocity, and the discharging time of the x-th stream, 
that of course depend on the thickness h of the brick. Because of physical 
considerations, we should expect that smaller is the diameter of the stream, 
i.e. less is the parameter x (or larger is a), smaller is the value V{x,h) of 
the velocity of the flow because more relevant is the effect of the friction 
forces and, consequently, larger is the discharging time T{x,h). It is worth 
to stress that the assumption of small cross section, that allowed us to write 
the velocity along the streams as in eq.(7), doesn’t affect our results form 
a qualitatively point of view, but has only been made in order to solve the 
integral in eq.(6) in a closed form. Substituting eq.(7) into eq.(6), we get the 
following expression for the effective velocity: 

'yefr(t) = :j- [ dxx^^~^^V{x,h)[H(t) — H(t— T{x,h)] (8) 

log c Jo 

In the limit of small cross sections, we can Taylor expand the x—dependence 
of V{x, h) and so write, to the lowest signihcant order, V{x, h) = V{h)x^ 
and T{x,h) = h/V{x,h) = {h/V{K)) x~^ = f{h)x~^ (with a coefficient 
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defined from a proper physical model or experimentally deduced, and T{h) = 
h/V{h)). Substituting such relations in eq.( 8 ), we then have 

Vesif) = j [H{t) - H{t - T{x, h)] (9) 

Such integral can be carried out noting that the integrand is not null only 
in the time interval 0 < t < T{x,h), obtaining: 

V'(h) 

Vesit) = n , t for f>0 (10) 

h loge 

in which we indicated 7 = l + (2 — Zl)//5 and V'{h) = V{hY~'^ W. Eq. (10) 
proves that the time evolution of the overall flux is ruled by a power-law, 
whose order is related to the anomalous dimension D of the cross section of 
the brick, and to a coefficient /9, characteristic of the model used to describe 
the motion of the fluid along the streams. 

Let us introduce now a physical model into the general arguments outlined 
so far. As anticipated, we model the energy losses in the flow due to the 
friction forces, by using the Darcy-Weisbach equation, according to which 
such losses, after a length L, per unit mass of the fluid , are equal to: 

AEioss = (11) 

2g d 

where d is the “hydraulic diameter” of the stream’s section, and K is a 
dimensionless coefficient called the Darcy friction factor. The value of this 
factor can be deduced from the Moody’s diagram and, in the assumption of 
laminar flow (that is the limit we are working on), it is equal to K = 64/i?e, 
being Re = — the Reynolds number of the flow. The Bernoulli equation, 
by which we describe the motion of the fluid, particularized for the problem 
at hand (see fig.3), reads as 


-+ 64 


u \ {h — y) 


pv d J 2g d 


= s(l‘- p) 


( 12 ) 


having indicated by y the distance of the surface of the column of fluid, in 
the generic stream, from the top of the brick, by p and u respectively the 
density and the dynamic viscosity of the fluid, and hy g = 9.81m/s^ the 
standard gravitational acceleration. Solving this equation for the velocity v 
of the fluid, we obtain: 
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V 
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v{d,y) 


-32 


V 

pgd? 


{h-y) + 


pgd^ 


{h - yY + 2g{h - y) 


(13) 


so that the time instant at which the snrface of the colnmn of water is at 
the distance y from the top, is formally given by: 



(14) 



Fignre 3: Discharging of the viscons flnid throngh the hollowed brick, dehned 
in the text as a simplifying model for the Lichtenberg hgnre reported in hg.l. 


Becanse of the assnmption of small area of the streams, and so small 
hydranlic diameters, the eq.(13) can be approximated to the lowest order as: 


, , 9V<i" ,g‘pP 

gyy) ~ 


(15) 


32z/ 32z/ 

being ^ a coefficient relating the hydranlic diameter and the characteristic 
length £ of the streams. Remembering that snch length scales as i{x) = iox, 
we can identify in eq.(15) the parameters V(h) = = 2, previonsly 

defined. Inserting in eq.(lO), the Hansdorff fractal dimension D = 1.8928 for 
the sqnare Sierpinski carpet, and D = 1.5849 for the triangnlar Sierpinski 
carpet, together with the valne (3 = 2 calcnlated above, we obtain that, for 
snch geometries of the cross section of the brick, the flnxes are proportional 
respectively to: 4>sq{t) ~ 0Tr(^) ~ being the former rela¬ 

tive to the sqnare carpet and the latter to the triangnlar carpet. 
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In order to verify these results, and so the validity of eq.(7), we have calcu¬ 
lated numerically and directly from the eqs.(4), (13) and (14), the flux as a 
function of the time, by using the values e = 3, iV = 8 that dehne the square 
carpet and e = 2, N = 3 that dehne the triangular one. We assumed the 
side = 1 cm for the main stream, and a thickness h = 1 m of the brick. 
Furthermore we considered the density p = 1000 Kg/m^, and the viscosity 
V = 0.001 Pa • s for the huid. Solving numerically the system composed by 
eqs.(13) and (14), we computed the velocity of the how as a function of time, 
for fractal orders ranging between 0 < a < 4, with a step Aa = 0.2. The 
result is shown in hg.4, and reveals the goodness of the approximation we 
made, in writing v{x,t) in terms of step functions. Considering now the 



time (s) 


Figure 4: Velocity as a function of time, for pipes from the order a = 0 to 
a = 4, with a step Aa = 0.2 

wider range from a = 0 to a = 6, and a step Aa = 0.1, we obtained the 
time evolution of the overall hux from the brick, shown in log-log scale in 
hgs.5(a) and 5(b) respectively for the square and triangular carpet. Apart 
for a deviation in the long time domain, due to the hnite number of fractal 
orders used in the calculation, the time dependence is clearly a power-law. 
By htting these numerical results, disregarding the data relative to the long 
time domain, we get that the hux, for the square and the triangular carpet 
respectively, is proportional to (/>sq(i) ~ t~^-057 0Tr(^) ~ results 

that are very close to the ones calculated from eq.(ll). 
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(a) Square Sierpiknsi carpet 
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(b) Triangular Sierpiknsi carpet 

Figure 5: Overall flux as a function of time 
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4 The rise of fractional operators in the linear 
regime: the seeping problem 

In the previous section we dealt with the problem of percolation of a fluid 
through a porous medium, that has been modelled by the fractal hollowed 
brick. We showed that, as the flow takes place on every scale of the fractal, 
then the overall out coming flux evolves in time according to a power-law, 
whose order keeps memory of the scaling properties of the underlying fractal 
geometry, as well as of the physical model by which the motion of the fluid 
has been described. Basing on this remarkable result, our goal now is to give 
evidence of the occurrence of fractional operators in the equations ruling the 
evolution of such flow, once a pressure history is applied on the medium. We 
will show in particular that it happens once the response of the system to 
the forcing action is linearised around a hxed value, condition that allow us 
to take advantage of the Boltzmann superposition principle. This is exactly 
the same hypotesis on which the well known theory of linear viscoelasticity 
is based. As long as the value of the load F applied to a specimen of vis¬ 
coelastic material is sufficiently small in fact, then the non linear dependence 
of the elongation u on the load itself, found by Nutting in 1921 in his famous 
experiment [41], and equal in general to 


u = arF^ (16) 

can be linearised around the value F = 0. In eq.(16) m and n are parameters 
characteristic of the material, that vary with the temperature, while the 
coefficient a depends on the type of experiment performed. Again, in this 
assumption the Boltzmann superposition principle is now valid so that, as we 
consider the specimen subjected to a load history, the time evolution of the 
elongation is ruled by a linear fractional differential equation. By following 
the same arguments here described for the viscoelasticity, we will linearise 
the dependence of the flux on the pressure applied to the brick, obtaining a 
relation between these two quantities, of non integer order. 

Since it is not signihcant in this context, we disregard in what follows 
the effect of the gravity and of any other volume forces acting on the fluid, 
assuming it subjected only to the friction forces and to the net pressure acting 
on the two sides of the brick (that we assume empty for t < 0). At the time 
instant t = 0, we imagine to connect it to a reservoir, in which the fluid 
is at rest and whose pressure is assumed to be constant in time. Our hrst 
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aim is then to derive an expression for the time evolution of the flux seeping 
through the brick, as a function of the net pressure p applied to the brick. 
Again, in the assumption of sufficiently small area of the streams, the friction 
forces are predominant respect to the pressure ones and so the velocity in 
every streams can be assumed constant. In this hypothesis we can then write: 


v'{x,t) = V{x,p)H{t — T{x,p)) (17) 

where V{x,p) and T{x,p) are, respectively, the velocity of the fluid at the 
lower side of the brick, and the time necessary for for the fluid itself to cover 
the distance h, equal to the thickness of the brick. Inserting eq.(17) into the 
eq.(6) for the effective velocity Vesit), we get 


^eff (^) 


1 

- - / dxx'^^~^^V(x,p)H(t— T(x,p)) 

loge Jq 


(18) 


Describing again the motion of the fluid by using the Bernoulli equation, 
and the energy losses by the Darcy-Weisbach law, because of the hypothesis 
of small cross section of the streams, we can Taylor expand the expressions 
for the velocity and the discharging time, obtaining, to the lowest signihcant 
order, V{x,p) = V{p)x^ and T{x,p) = h/V{x,p) = h/V{p)x~^ = T{p)x~^, 

with V{p) = Vop = p, f{p) = fo/p = (^ 1 ^) I and again with 

(5 = 2. Since the unit step function in eq.(18) is not null for t > T{x,p) = 

/ N -1//3 

T{p)x~^, condition that is satisfied for x > ( i the integral can be 

solved in a closed form, giving 


t'efr(t) 


vjp) V" 

log 6/^7 \T{p)) 


for t > T{p) 


(19) 


where, as before, ■y = 1 + {2 — D)/(3. We obtained again a power-law ruling 
the evolution of the flux, whose order is exactly the same as the one obtained 
in the previous section. This remarkable result conhrms that the order of the 
power law depends only on the anomalous dimension of the underlying fractal 
geometry and on the model used to describe the physics, and is independent 
on the particular situation at hand. As done in the previous section, we 
verihed the validity of eq.(17), by calculating numerically and directly form 
the Bernoulli equation, the flux as a function of time. We considered the 
range of fractal orders between a = 0 and a = 10, with a step Aa = 0.05. We 
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isolated the time dependent part of such solution, that is reported in £g.6(a) 
and 6(b) in log-log scale. Disregarding the short and long time domains, in 
which a deviation from a power-law is evident and due to the hnite number 
of the fractal orders used in the calculation, we htted the data, obtaining the 
time dependences for the flux ^sq ~ ^ ^- 1 . 2075 ^ respectively 

for the square and triangular Sierpinksi carpets, that again are very close to 
the one obtained from eq.(19). 



(a) Square Sierpiknsi carpet 



(b) Triangular Sierpiknsi carpet 

Figure 6: Overall flux as a function of time 

So far we have calculated the response of the system (in terms of the 
overall flux of the viscous fluid) to a step-like excitation. In section 2, such 
excitation was given by the gravitational forces that, switched on at the time 
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instant t = 0, enforced the fluid to flow through the medium. In this section 
instead, we dealt with the case in which the flow is driven by a constant net 
pressure applied to the medium, again for t > 0. We found in both cases, 
a power-law for such response, whose order is related to the dimension of 
the fractal shape by which we modelled the pattern of the streams along 
which the fluid flows through the medium. Basing on these results, the last 
fundamental step we want to perform, and that represent the principal aim 
of the paper, is to give evidence of the emergence of fractional differential 
equations, in ruling the response of the system, once it is subjected to a 
load history. Because of the linearity of fractional operators, we prove the 
validity of this statement, once the non linear dependence of the response on 
the amplitude of the forcing action, is properly linearized. In particular, in 
the present context, we prove the emergence of fractional operators in the 
constitutive law relating the overall flux of the fluid and the net pressure 
history applied to the brick, once the eq.(19), that represent the response 
of the system to a step-like forcing action, is linearized around a hxed value 
p = Po- To this aim, let us assume such pressure history varying in a small 
region of pressure values around Pq. This assumption allows us to linearize 
eq.(19) in the neighbourhood of p = pq. Disregarding the constant pressure 
contribute, and focusing on the effect due to Ap(t) = pit) —po, by exploiting 
the Boltzmann superposition principle, we can write the effective velocity in 
the form: 


Cl 


Ap{t) + C 2 


r‘t-T(po) 


Ap{t) 
{t - r)7 


dr 


for t > T(po) (20) 


where we dehned the coefficients 


Cl 

C2 


% 

h log e 

(7-1)ToW 


( 21 ) 


In the eq.(20) for the effective velocity, an integral term appears, that looks 
like a fractional operator. It is easy to show that integrating it by parts, 
and assuming the initial conditions Ap{t = 0) = 0 and Ap{t = 0) = 0, the 
constitutive law relating the effective velocity (and so the overall flux) and 
the pressure history applied to the brick, is a differential equation of not 
integer order that reads, again for t > T{p), as: 
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Weff(t) = CiAp {t - T{po)) + (t - T{po)) 


-r(2-7) 


Co 


7 - 1 


D 


f{po) 


Ap (t-f(po)) (22) 


where (t) is the Caputo’s dehnition of fractional derivative of 

order 7. The eq.(22) is the fundamental result we were looking for. It proves 
that, in the hypothesis of linear regime, a fractional differential equation nat¬ 
urally comes up in ruling the evolution of a natural phenomenon or process, 
taking place on an underlying fractal geometry, whose order is related to the 
anomalous dimension of the latter. We derived such result by studying the 
physics of the flow of a viscous fluid through a fractal shape. Anyway the 
considerations made, that lead us to such remarkable result, are general and 
so hold true whichever is the physical problem at hand. 

This result gives us a deep insight into the world we live in. We know 
that differential operators of integer orders, as well as integer order ordinary 
or partial differential equations, are peculiar of smooth geometries, on which 
all the physical theories developed up to now, from the newton laws of me¬ 
chanics, to the Einstein theory of gravitation, or quantum mechanics, are 
based. On the other hand, according to the results obtained in this paper, 
the success of fractional calculus in modelling physical phenomena, reveals 
an underlying fractal nature of reality. Because of that, we are conhdent that 
in the next future, the fractal geometry will have a leading rule in gaining a 
better and more complete understanding of physics and of the others natural 
sciences, hopefully clarifying some of the issues by which are still affected. 


5 Conclusions 


In this paper, we showed that a relation between fractional calculus and frac¬ 
tal geometry exists, which is intimately related to the physical origins of the 
power-law long memory and hereditary properties observed in many natural 
phenomena, and that are the characteristic feature of fractional operators. 
With the aid of the relevant example of a viscous fluid seeping through a 
porous medium, that we modelled as a fractal shaped hollowed brick, we 
showed that a power-law naturally comes out, ruling the evolution of a phys¬ 
ical phenomenon, every time it takes place on an underlying fractal geometry. 
Trace of the scaling properties of the underlying fractal geometry, exhibits in 
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the order of such power-law that, accordingly, we found to be strictly related 
to its anomalous dimension, and to a coefficient, characteristic of the model 
by which we describe the physics of the phenomenon at hand. 

Linearizing the non linear dependence of the response of the systems at hand 
to a proper forcing action we proved, by exploiting the Boltzmann superpo¬ 
sition principle, that a fractional differential equation naturally comes up in 
ruling its evolution. The order of such equation turned out to be related to 
the anomalous dimension of the geometry on which the phenomenon takes 
place. 
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